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Поставлена и решена начально-краевая задача системного исследования концентрационного и электриче- 
СКОГО ПОЛЕЙ В двухэлектродном электрохимическом устройст, ве. Рассмотрены вопросы сопряжения этих по- 
Лей. Получена электрическая схема замещения электрохимического уст. ройст, ва как объекта с распределен- 
ными параметрами, позволяющая вести строгие расчеты систем, содержащих эти ус тройс тва. 
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Введение. При строгих расчетах и исследовании динамики электрохимических устройств требу- 
ются их современные математические модели. Последние могут быть получены только в резуль- 
тате системного исследования в них процессов электро- и массопереноса в динамическом режи- 


и(0 








Рис. 1. Геометрия 
электрохимического устройства 


ме. Важным моментом этого исследования является замкнутая 
связь концентрационного и электрического полей с помощью 
процессов, протекающих на границах раздела «электрод - элек- 
тролит» — электродных процессов. Конечной целью системного 
исследования электрохимического устройства является получение 
его математической модели как элемента электрической цепи. 

Постановка задачи. Исследуется электрохимическое устройст- 
во, состоящее из двух одинаковых плоскопараллельных электро- 
дов. Через это устройство проходит ток Г (8, причем плотность 
тока по поверхностям электродов распределена равномерно. 
Площадь электродов $, расстояние между электродами # Лимити- 
рующей стадией кинетики электродных процессов является моле- 
кулярно-гиперболическая диффузия в электролите, протекающая 


с конечной скоростью И = \/Б/х, , где О -— коэффициент диффу- 
зии, т, — постоянная релаксации. Будем считать пространствен- 


но-временное концентрационное поле электролита с (х; 9 одно- 
мерным, причем координата х нормальна поверхностям электро- 
дов. Электродам соответствуют координаты х = дих= #(рис. 1). 


Относительно концентрационного поля электролита с (х; 8 ставится [1] следующая на- 
чально-краевая задача на отрезке [0; Й: 


дс 9?с 0?с 











тт, р: р 2! (1) 
О ОЕ 9х 
с(х; 0) = с, (2) 
О Ве 
5=(%:0) =0, (3) 
Е Ц 
5 (020) М =. (4) 
Об вы) 
(мг, (5) 


где М> 0 - кинетическая константа электродной реакции; с, - начальная концентрация элек- 


тролита. 
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Системное исследование полей. Задачу (1) - (5) удобно решать операторным методом Лапла- 


са. Пусть имеют место соответствия с(хзе |= (х;р), ЦКВ-Гр). 


Относительно изображения с(х; р) получаем следующую краевую задачу: 











ас(х;р) _ р(ч.р +1). с (&р+ 
2 И о О. (6) 
Вер) АМ И. (7) 
Хх 5 
ве Р) ур). (8) 
Хх $5 


Общее решение дифференциального уравнения (6) имеет следующую структуру: 
с(х;р) =6(х;р)+с,(х;р), (9) 
где (х; р) — общее решение соответствующего однородного дифференциального уравнения 


а26(х; О 
си р РР 2 (ж;р) =0, (10) 





Сн(х;р) - некоторое частное решение исходного неоднородного уравнения (9). 
Для решения (10) составим характеристическое уравнение 


<2_ Р.Р + _(, (11) 
р 


ай [2(<,р +1) 
ур — 


позволяют записать общее решение (10) в виде: 


без) = арен Ри + веру М, (12) 


где А (р), В (р) — произвольные постоянные интегрирования, подлежащие определению. 


Корни уравнения (11) 


Частное решение с.(х; р) уравнения (6) можно записать в виде: 
с.(х;р) = 6. (13) 
р 
Подстановка (12) и (13) в (9) приводит к такому результату: 


С(хур) = А(р)зн РР +В + в(русн, Р-Р +В, ©. (14) 
р р р 


Коэффициенты А (р) и В (р) находим из краевых условий (7) и (8). 
Для этого продифференцируем (14) пох: 


Чс(ж:р)_ [2@,2+1) [Р.Р +1) [Р(.2+0 
ро р | ож р х-+В(р)58 р *] (15) 


Из (15) прих = 0 имеем: 


4с(0;р) _ у [Р.Р +1) 
А Ар, (16) 
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Сопоставление (7) и (16) приводит к такому результату: 


А(р)=м— КР) (17) 


с Ра. 
р 


Подстановка (17) и х = в (15) дает: 


4с(вр) _ МИР), [Р.Р +1), О 2. +В ‚Е +1), 18) 
ах 5 р | р р ” 


Сопоставление (8) и (18) определяет коэффициент Бр). 


Е (+р+1) 6 
Грзв. РР © 
В(р)=М НЕЕ (19) 


5 Вер 4 
р р 2 























Общее решение с(х;:р), заданное выражением (14), с учетом (17) и (19) приобретает 


следующий вид: 
№ Г(р)5В Р(ер +1) (х :] 
: я у р 2 
с(х;р)=— + 


р 5 Эт ер +В й 
р р 2 
В любой электрохимической системе на границе раздела «электрод — электролит» сущест- 
вует скачок электрического потенциала, однозначно определяемый при принятых нами допуще- 
ниях значением концентрации электролита на граничной поверхности. Эта зависимость задается 
уравнениями Нернста, которые в исследуемом случае задаются следующими линейными соотно- 
шениями [1]: 











(20) 





А (В) = 9, + 9,6(0;#), (21) 
А* (Е) = 9, + 9,С(6;1), (22) 
где Л (ГР) - скачок потенциала на катоде; л*(РГ) - скачок потенциала на аноде; 


9 >0, 9, > 0 - параметры линейной аппроксимации уравнения Нернста. 


Применяя преобразования Лапласа к соотношениям (21) и (22), получаем следующие 
операторные зависимости: 


^-(р)= г +9.(0;5), (23) 


А (р)= 9,6). (24) 


Напряжение на электрохимическом устройстве (/ (#) рассчитывается с помощью 2-го зако- 
на Кирхгофа: 


А* (Е) + (г - А -И( 50, 


И 
где Г =—— - сопротивление столба электролита между пластинами. Здесь у. - удельная 


и 


электропроводность электрона. 
Из 2-го закона Кирхгофа получаем: 


ИВ = А -^ (+ Це. (25) 


((р) = (р) - (р) + Круг, (26) 
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Подставим (23) и (24) в (26). 
Ц(р) = 9,(с(6;р)- с(о;р)) + Г(р)г,. (27) 


Получим с(о;р) и с(е:р) из (20): 


зн. Р-Р +1 4 
МГ(р). 2—2 (28) 


5 ор в й" 
р р 2 


сн Р-Р +1) 6 
р 5 [РР + ГР,р-+1) с 
р р 2 


Из (27) с помощью (28), (29) получаем: 


2ма, зн, | РР +1) - 
(р) = Р2 Кр)+кЦр.. (30) 


Е "Веер +1) [ 
р р 2 


Синтез электрической схемы замещения. Полученное соотношение (30) лежит в основе син- 
теза электрической схемы замещения. 








с(о;р) = + 
р 




















Коэффициент при 1 (р) в первом слагаемом (30) является диффузионно-гиперболическим 


импедансом: 
т ы 5 
2(р)= 


52 + бас +1 
р В. 9 


Обратная ему величина - проводимость — выражается так: 


(р) -_5 ен |еР+ 
20, р В: 0 


Разложим [2] гиперболическую функцию, входящую в (31), в следующий ряд: 


а [Ра,.р+1) = 2 р (р(«,р +1) | — у 1 , 
р 2 р(х р- 1) р ЩЖ к Р(‹,р +1) 2 2 
мм ИИ +К 
Гр) 4)т 


Подставляя последнее разложение в (31) и опуская промежуточные выкладки, получаем: 
25 › 25 

















(31) 














5 ЗМР 09М 
У(р) = +29 В 32 
(2) 2№9, 2 2. ‚ 1_ 42 ВЯ 
Р ее с. 


Проводимость У(р), определенная формулой (32), отвечает условиям [3] физической 
реализуемости пассивными электрическими элементами. 

Покажем, что каждый член ряда в формуле (32) можно моделировать электрической схе- 
мой (рис. 2). 

Импеданс схемы, приведенной на рис. 2, имеет вид: 

2 (р)- РС, (+) + р + СГ) + г, 
К 2 . 
РС, (рЬ, + г.) 
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Рис. 2. Электрическая схема ветви замещения 


Соответственно, проводимость ветви определяется выражением: 
2. 
РС, (рЬ, + г.) 








у = : 33 
«(р РС, (О +) + р + С гг) + к. 9 
Приведем выражение (33) к виду (32): 
м 
(1) (2) (1) (2) 
У, (р) = АВЕ, ВЕР (34) 





(и +г,0)) г. @® 
Бо”): ВСЕ) 
Сопоставляя соответствующие коэффициенты в (32) и (34), получаем систему уравнений 
для определения параметров ветви электрической схемы замещения: 
1 Е. 
И Е 





р’+р 











в” 28$ 
С, (г, + г.) 609, Мт, ' 
НВ. 
С, (и + и) т, ' 
в” 4К??) 





ВС +к о) = т, 
Решение последней системы уравнений приводит к следующим результатам: 
$4 х_ 9/4 о 22,9 МО 2К°л’ 09. М 
С, = ар = и = - ‚ Ш= - : 
29. Мк) 251 5( 54 
Учитывая структуру выражений (32) и (30), получаем полную электрическую схему заме- 
щения электрохимического устройства (рис. 3). 











С, 






С 


5) и 





о п) 
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Рис. 3. Электрическая схема замещения электрохимического устройства 
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В этой схеме сопротивление г, , моделирующее омические потери процесса поляризации 


электродов, вычисляется по формуле: 
При исходных данных 
5 = 10м7"; ( = 0,01м; д, = 0,025 В/(кмоль/м?); № = 5,45 (кмоль/м3) мА”; 
т2 = 1108 с; у. = 35 Ом**м"; р = 1,65'1078 мс! 
получаем следующие значения параметров электрической схемы замещения: 
Ь = 2,857 Ом; г, = 13,625 Ом; 
С, = 0,225 Ф; п = 6,813 Ом; п® = 4,438-10-8 Ом; [1 = 4,438*10°° мкГн; 

С = 0,0563 Ф; 5 = 6,813 Ом; 2) = 17,752:108 Ом; [> = 17,752:108 мкГн; 

С; = 0,025 Ф; в? = 6,813 Ом; 3 = 39,942:10° Ом; (3 = 39,942:10°8 мкгГн. 
Заключение. Полученная математическая модель и электрическая схема замещения электрохи- 
мического устройства могут быть использованы для инженерных и научных расчетов режимов 
работы этих устройств. Особенно эффективным использование этой модели может быть при изу- 
чении релаксационных и ключевых (быстропеременных) режимов коммутаций. 
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